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Abstract 



We revisit the strong coupling limit of the Schwinger model on a lattice using staggered 
fermions and the hamiltonian approach to lattice gauge theories. Although staggered fermions 
have no continuous chiral symmetry, they posses a discrete axial invariance which forbids 
fermion mass and which must be broken in order for the lattice Schwinger model to exhibit 
the features of the spectrum of the continuum theory. We show that this discrete symmetry 
is indeed broken spontaneously in the strong coupling limit. Expanding around a gauge 
invariant ground state and carefully considering the normal ordering of the charge operator, 
we derive an improved strong coupling expansion and compute the masses of the low- lying 
bosonic excitations as well as the chiral condensate of the model. We find very good agreement 
between our lattice calculations and known continuum values for these quantities already in 
the fourth order of strong coupling perturbation theory. We also find the exact ground state of 
the antiferromagnetic Ising spin chain with long range Coulomb interaction, which determines 
the nature of the ground state in the strong coupling limit. 

lr This work is supported in part by the Natural Sciences and Engineering Research Council of Canada, the 
Istituto Nazionale di Fisica Nucleare and NATO CRG 930954. 



1 Introduction 



Anomalies occur when all of the symmetries of a classical field theory cannot be realized 
simultaneously in the corresponding quantum field theory . The classic example is the Adler- 
Bell-Jackiw anomaly [|TJ in a vector-like gauge theory where gauge invariance and global axial 
symmetries are incompatible. This is a manifestation of the fact that it is impossible to find 
a regular ization of the ultraviolet divergences preserving both axial and gauge invariance. If 
gauge invariance is preserved, the observable manifestation of the anomaly is the failure of 
the conservation laws for axial currents and the resulting absence of their consequences in the 
spectrum of the quantized theory . 

An interesting issue which arises in the lattice regularization of gauge theories is how the 
lattice theory produces the effects of the axial anomaly. Lattice field theory is manifestly gauge 
invariant by construction and therefore, to produce continuum theories, they must find a way 
to violate axial current conservation. Normally, in a lattice gauge theories axial anomalies 
are either cancelled by fermion doubling or else the lattice regularization breaks the axial 
symmetry explicitly ||. The lattice Schwinger model, in the Hamiltonian (latticized space, 
continuum time) formalism 0] with one species of staggered fermions|3]] ||, represents the 
unique example where neither of these occur. The e ffects of the anomaly are not cancelled by 
doubling since the continuum limit of a 1+1-dimensional staggered fermion produces exactly 
one Dirac fermion which is the matter content of the model (see Appendix A). Furthermore, 
even though the continuous axial symmetry is broken explicitly by the staggered fermions, a 
discrete axial symmetry remains. It corresponds to the continuum transformation 

ip(x) — >j 5 i/)(x), tJ)(x) — > -^(x)7 5 , (1.1) 

and appears in the lattice theory as a translation by one site. Since the mass operator ipip is 
odd under ( |1.1|) , if this symmetry is unbroken then < tpip >= 0. However, in the continuum 
theory it is known that M, 

#Mx)) = -^ . (1.2) 



where 7 = 0.577... is the Euler constant. The expectation value in eq. (|1.2|) is the chiral 
condensate. This non-zero value of the chiral condensate is one of the main effects of the 
anomaly. In order for the lattice theory to produce this and other features of the continuum 
theory, the discrete axial symmetry ( |1 . 1|) must be spontaneously broken. 

The continuum Schwinger model J7|, BJ is the prototypical example of a solvable model 
where the anomaly occurs. It is 1+1-dimensional electrodynamics with a single charged 
massless Dirac spinor field. The action is 

S = J d 2 x@(ij^ + 7m A")V - ^F^) (1.3) 



2 It is worth to stress that the sign of the chiral condensate in eq.(1.2) is just a convention. The chiral 
condensate has always the opposite sign of the fermion mass m that appears in the Lagrangian of the theory: 
it is the order param eter for a spontaneous symmetry breaking of the chiral symmetry. The chiral condensate 
is the analogous of the magnetization for a spin model and the mass m plays the role of the magnetic field. 
So even if one is studying a massless model, in order to dec ide the sign of the chiral condensate, one has to 
consider a small mass term m that then one sends to zero. The minus sign which appears in eq.([L^) is due 
to the fact that one considers the m — > + limit in the Schwinger model Lagrangian. 
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It is invariant under gauge transformations, 



il>(x) — * e ie{x) ij(x)^\x) — > ^(x)e~ ie{x) (1.4) 

^ ► Afj, + <9^(x) (1.5) 
and, formally, under the axial phase rotation, 

tp(x) — > e^ a tlj(x),^(x) — ► i)\x)e~^ a (1.6) 

(where 7 5 = i / y 0/ ~f 1 ). 

At the classical level the above symmetries lead to conservation laws for the vector and 
axial currents, 

f{x) = ^(rr) 7 ^(x), jg{x) = ^(xhVW, (1.7) 

respectively. At the quantum level both currents cannot be simultaneously conserved. If 
the regularization is gauge invariant, so that the vector current is conserved, then the axial 
current obeys the anomaly equation |J 

= ^F, u {x). (1.8) 

Moreover, the correlation functions of the model do not exhibit axial symmetry. 

There are only neutral particles in the spectrum of the Schwinger model; notably, the 
bound state boson created by the axial current operator, 

d»<f>(x) = y/Zjg(x) (L9) 

appears in the spectrum as a free pseudoscalar with mass 
In terms of the boson field, the action is 

S = J d 2 x{\d^(j> + ^4> 2 ) (1.10) 

and the correlation functions for the vector and axial currents are obtained from the correlation 
functions of derivatives of the bose fields. The coupling constant e c has the dimension of a 
mass and the model is super-renormalizable. 

Since the Schwinger model is solvable it has often been used as a field theory where 
methods of lattice gauge theory, particularly the strong coupling expansion can be tested and 



compared with exact results of the continuum model ||10||[pT|. 

In this Paper, we shall revisit the lattice quantization of the Schwinger model. Our main 
purpose is to demonstrate two results. One is to show how the effects of the anomaly appear 
through spontaneous symmetry breaking. The other is to demonstrate that a strong coupling 
expansion around a gauge, parity and charge conjugation invariant lattice ground state can 
produce results in good agreement with the continuum. Our results turn out to be quite 
different from those of previous strong coupling computations |[L0||||11|| . The difference arises 
from our careful treatment of the normal ordering of the charge density operator, which must 
be defined so as to be compatible with all of the discrete symmetries of the theory. The 
necessity of careful normal ordering of the charge operator was pointed out in |T2]] . Our 
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vacuum state, vacuum energy and energies of elementary excitations are different from those 
in pi ITT], [Tj| and, as we shall show, they give a more accurate extrapolation to the continuum 



limit. 

The Hamiltonian and gauge constraint of the continuum Schwinger model are 

H = J dx (^E 2 (x) + %l)\x)a{id x + eA(x))^(x)j (1.11) 

d x E(x) + i)\x)i)(x) ~ (1.12) 
A lattice Hamiltonian and constraint which reduce to these in the continuum limit are: 

Hs = ^T. E l ~ £EW4 + ie^ - V4e""-^i) (1-13) 

E x -E x ^ + 4^ x -^~0 , (1.14) 

where the fermion fields are defined on the sites, x = — Vj - y + 1, gauge and the electric 
fields, A x and E x , on the links [x; x + 1], N is an even integer and, when N is finite, we use 
periodic boundary conditions. When N is finite, the continuum limit is the Schwinger model 
on a circle [13|, ||. The coefficient t of the hopping term in ( 1.13Q plays the role of the lattice 



light speed. In the naive continuum limit, = e c and t = 1. However, we shall keep 6l and 
t as parameters which can be adjusted to fit the lattice values of quantities such as the mass 
gap and the chiral condensate to those which are known in the continuum. 

The non vanishing (anti-)commutators for the lattice variables are 

[A x ,E y ] = i5 xy , {i>x,i>l) = $xy (1-15) 

The Hamiltonian and gauge constraint exhibit the discrete symmetries 

• Parity P: 

A x — -+ -A.*-!, E x — > ip x — > (-l) x V-x, 4 — > (-l)^-x (1.16) 

• Discrete axial symmetry T: 

A x — ► A x+1 , E x — ► E x+1 , i) x — ► ip x+1 , ipl — ► ipl +1 (1.17) 

• Charge conjugation C: 

A x — > -A x+1 , E x — ► -E x+1 , ip x — > ipl +1 , ipl — ► ip x+1 (1.18) 

When A x = 0, the spectrum of the hopping Hamiltonian is e(pa) = tsmpa(p G [0, — ]). 
In the low energy limit it resembles a massless relativistic spectrum for excitations near 
the Fermi level. The two intersections of the energy band with the Fermi level provide 
the continuum right and left moving massless fermions. The electron ground state is 
invariant under charge conjugation only when the Fermi level is at e(j)p) = (i.e. when 
exactly half of the fermion states are filled and J2 X < p(%) >— 0). In the remainder of 
this paper we consider only this case of half-filling. 
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One has two possile ways to put the spinors ip 



in 

^2 



on the lattice in the staggered 



fermion formalism. One can put upper components on even sites and lower components on 
odd sites or viceversa. The mass operator, which reduces to i/ji/j in the continuum limit is the 
staggered charge density 

(-IY 



M{x) = 

for the first choice of staggered fermions or 



tl>x - 1/2 



M{x) 



- 1/2 



1.19) 



1.20) 



in the other case. 

These operators ( |1.19| ) ( P^2H| ) are scalars under parity and are even under charge conjuga- 
tion but change sign under a discrete axial transformation. 

The lattice Schwinger model is equivalent to a one dimensional quantum Coulomb gas on 
the lattice. To see this one can fix the gauge, A x = A (Coulomb gauge). Eliminating the non- 
constant electric field and using the gauge constraint, one obtains the effective Hamiltonian 



He 



+ 



H u + H p = 



it 



4 £ 2 + ^L^p{x)V{x-y)p{y) 



2N 



x,y 



+ 



where the charge density is 
and the potential 



p(x) 



V(x - y) 



N 



N-l 
n=l 



i2irn(x—y)/N 



4 sin 2 £2 



1.21) 



(1.22) 



(1.23) 



is the Fourier transform of the inverse laplacian on the lattice for non zero momentum. The 
constant electric field is normalized so that [A, E] = i . The constant modes of the gauge field 
decouple in the thermodynamic limit N — ► oo. The gauge fixed Hamiltonian Ql.21 ) can also 
be written as a quantum spin model. Consider the Jordan- Wigner transformation 



U[2iS 3 (y)]S-(x) 

y<x 

Then, the Hamiltonian is 

tcos(A) 



iPi=Yl [-2iS 3 (y)]S + (x) , VdV. = S 3 (x) + 1/2 (1.24) 

y<x 



H 



(S 1 (x)S 1 (x + 1) + S 2 (x)S 2 (x + 1)) 4 



e?a 



fJ2S 3 (x)V(x-y)S 3 (y)+elE 2 /2N 

x,y 

(1.25) 



For the moment, we ignore the constant modes E and A. The first term in the Hamiltonian 
is the quantum x-y model which has a disordered ground state (of course, exact solution of 
this model is obtained by the inverse Jordan- Wigner transformation [|14], |I5|). The second is 
a long-ranged Ising interaction. We will prove (see also ref. ||12|| ) in the next section that, 
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in spite of the low dimensionality of this system, the latter term has a Neel ordered ground 
state. This is a result of the infinite range of the Coulomb interaction. When both terms 
are present in the Hamiltonian, they compete, the first one favoring disorder and the second 
one favoring order. Later, when we extrapolate our strong coupling results to the continuum 
limit, we shall assume that the ordered ground state persists for all positive values of the 
dimensionless constant e 2 L a 2 /t. The strong coupling limit is where this constant is large and 
the ground state has Neel order. The continuum limit of the Schwinger model is found where 
this constant is small. We shall assume that the Neel order persists in that limit. f\ 

There are two interesting questions which we shall not pursue in this paper, but will 
consider later. One is the possibility of finding exact solutions of ( 1.25|) with E and A set 
to zero using the Bethe ansatz. This would effectively give an exact solution of the lattice 
Schwinger model in infinite volume. The other is the solution of the model ( p.. 25 ) when the 
topological modes, E and A, are included but the Coulomb interaction is ignored. The solution 
of this model would give some insight into how the anomaly arises in the weak coupling limit 
of the lattice theory. 

In section 2 we set up the strong coupling expansion around the gauge invariant ground 
state of the unperturbed Hamiltonian H u . Then, we compute to the fourth order in powers 
of e = „ t 2 , the corrections to the energies of the ground state and of the low-lying bosonic 

excitations (at zero momentum). Subtracting the energy of the ground state from those of 
the low-lying excitations defines the lattice meson masses. In section 3 we compute the chiral 
condensate to the same order in e. Section 4 is devoted to the comparison of the lattice results 
with the continuum ones; there we shall see that, having used a gauge invariant ground state, 
the improved strong coupling expansion leads to a by far more accurate evaluation of the 



physical parameters then the strong coupling expansions given in QTOJ, [IT], [Lq| . Section 5 
shortly comments on our results. 

In the appendices we provide some technical details. 



2 The low-lying excitation spectrum 

We shall examine the Schwinger model in the strong coupling limit. To solve the lattice 
Schwinger model in the strong coupling expansion, it is necessary to find states which are 
annihilated by the generator of gauge transformations 

E x -E x ^ + i)li) x -- (2.1) 

and, at the same time, are eigenstates of the unperturbed Hamiltonian 

Hu= 6 fEE 2 x . (2.2) 

^ X 

We shall quantize the gauge fields in the functional Schrodinger picture where the wave- 
functions are functions of the variables A x and the electric field operators are defined as 

E x = \— - . (2.3) 

% 8 At ^ ' 



3 If it did not persist, but there was a second order phase transition at some critical value of the constant 
to a disordered state, the continuum limit would be found at that phase transition. 
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In a previous publication [|12j it was proved that there were two gauge invariant ground states 
of the system. In fact, when the number, N, of lattice sites is even the two states 



\tfj> = [J 4|0 > e^*^^ (2.4) 

x=even 

\ x > = n '^.|0>e i ^^^ (2.5) 

x=odd 

are degenerate gauge invariant ground states of H u . To see this, note that then H u acts 
on gauge invariant states, it is effectively the sum of two operators, the zero mode energy 
e 2 L E 2 /2N where E = Y,x Ex/ VN and the Coulomb energy 

^ £p(*)nr - y)p(y) = 5^E „(-*) 1^*) ■ ( 2 -6) 

z x,y Ziv fc^o 4sm 2 

where we have used the Fourier transform 

V4V>. -^ = ^E e*"p(p) ,P=j,^ {-f , f - 1} (2-7) 
Furthermore, gauge invariant states have vanishing total charge, 

YM>* - \) = ( 2 - 8 ) 

The Coulomb energy satisfies the bound 

e j^L y p < k )_L^ p <k) >^y p (-k)-p(k) = 

= e ^U^-b 2 = ela/32 (2.9) 
o x 2 

(since (ipli/} x — \) 2 = 4). The states |^ > and |x > are eigenstates of the charge density with 
eigenvalues p{k) = ±^5k, ir. These charge densities saturate the lower bound for the Coulomb 
energy. Furthermore, E\ip >= = E\\ > so \ip > and |% > minimize the energy of the zero 
mode. Therefore, they are degenerate ground states of H u . 

The ground states are invariant under parity and charge conjugation but they both break 
chiral symmetry, i.e. the symmetry under translations by one site. The order parameter of 
the chiral symmetry, the mass operator M = jr- J2x M(x), is diagonal in these states and has 
eigenvalue ~. In order to have a chiral condensate of negative sign as in the continuum, one 
has to consider the mass operator M'(x) = —M{x). It is interesting that the ground states 
are degenerate and are not invariant under chiral symmetry, even in finite volume and in one 



spatial dimension ; it is possible due to the long-ranged nature of the Coulomb interaction ||17||. 
The chiral symmetry is spontaneously broken in the infinite volume limit. It is essential that 
the number of sites of the lattice, N, is even. 

Of course once perturbations from the hopping term in the fermion Hamiltonian are taken 
into account, when the volume is finite, there are high order terms (of order N) in the strong 



4 One has to consider M(x) given by eq.(1.19) when applying it to \ip > and M(x) given by eq.( 1.20| ) when 
applying it to \\ >■ 
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coupling perturbation theory which mix the ground states that are constructed by perturbing 
\ip > and \x >■ The mixing amplitude is 

1 iV 2 \t\ N 4 JV ~ 1 AT' 

Since the amplitude ( |2.10| ) is negative and real , whatever the diagonal elements are, the state 
with lower energy is the state 

^(|V>>+|X>) (2.H) 

This is an eigenstate of the discrete axial transformation, with eigenvalue 1. Thus the true 
ground state is chirally symmetric. 

In order to avoid this restoration of the symmetry, it is necessary to take the infinite 
N limit so that the two ground states are not mixed to any finite order of perturbation 
theory. Thus, in the perturbative expansion one only has to consider diagonal matrix elements 
and, consequently, perturbation theory for non-degenerate-states. We shall then choose as 
unperturbed ground state \ip >. The results concerning the energies would have been the 
same had we chosen |% >. 

Excitations are created by the operators 



^yEL.^ (2.i2) 



and have energy — j\ greater than the ground state energy. A single electron, in order 
to have a gauge invariant state, must be attached to a line of electric flux which goes from 
the electron to infinity. This is not allowed at all in a finite volume with periodic boundary 
conditions and on the infinite lattice, the string of electric flux would have infinite energy. 

e 2 

Thus electric charges are confined and the string tension is given by the electric charge -y. 

Let us now start the perturbative computation of the excitations mass gap. We shall work 
in the Coulomb gauge V A = 0, so that A is x-independent and the ground states are 

> = n v£io> ( 2 - 13 ) 

x=even 

\X> = II V£|0>. (2.14) 

x=odd 

The Coulomb potential reads 

v{*-v) = ± E d?!^" 9 ' = w ( * - ^ - > - y i + m ^ - ^ (2 ' 15 > 

The Schwinger Hamiltonian, rescaled by the factor e|a, may then be written as 

H = H + eH h = -^(Hc + H p ) (2.16) 

e L a 

with 

#o = E - i -^ 1 WMy) (2.i7) 

x>y ZJV Z 



7 



Hh = -*£(V4W<d - ^ x e~ iA ^ x+l ) = -i{H R - H L ) (2.18) 



X 



and 

t 



£ =2i^ < 2 - 19 > 



To the fourth order in e the perturbative expansion for the energy is given byQ 

^ = Ef + ^ E f + ^ 



E, = Ef + e*Ef + e* E f (2.20) 



where 



Ef =< mo\i> > (2.21) 
4 2) =< ^\ H I tW^T H ^ > (2-22) 



^t-^ H h \ib >< iblHt — ^J- 



< ml- W ± —H h \^ >< m l, 0) * *W > (2-23) 



and 1 — n^, is a projection operator which projects a generic state \s > on \if> > 

(l-II^)|s>=|V> • (2.24) 

The phases e lA in the hopping Hamiltonian are irrelevant until we reach the N th order in 
the perturbative expansion. As a matter of fact, only at this order the transverse electric field 
becomes important. By direct evaluation one gets 

Ef = | (2.25) 

Ef = -AN (2.26) 

Ef = 192N (2.27) 

E^ = — — 4iVe 2 + 192iVe 4 . (2.28) 
32 



5 We observe that for the ground state perturbative expansion the relative minus sign in the hopping 
Hamiltonian is irrelevant whereas it will play an important role in the calculation of the excited state energy. 
This is easily seen if one notices that at the second order in e one gets 

Ef = <^\{H L -H R ) ^ (H r -H l )\tP>= 

- H 

= <4,\H L ^f^ H R \^> + <il>\H R -^!- H L \^> 

E,u -H E, b -Ho 
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Next we compute the energies of the low-lying excitations to the same order in the per- 
turbative expansion. The lowest -lying excitations are degenerate and are created by the 
operators 



H L (x) = ipte iA iP x+ i (2.29) 

c+l ( 



H R (x) = ^ x+l e~ lA ^ x (2.30) 



which move a particle one lattice spacing to the left and right respectively. When they act 
on the vacuum \ip > they create the two states 

\x,R> = H R (x)\ij >= iPl +1 e- iA iP x \iP > (2.31) 
\x,L> ee H L {x)\^>=^ x e iA ^ x+l \^> (2.32) 

which are degenerate in energy (see eq.( |7.2| ) in appendix B). In order to resolve the degeneracy 
of these two states , one must find the first nontrivial order of the perturbative expansion in 
which the two states are mixed; this turns out to be the second order. It is easy to show that 
the second order of perturbative expansion is diagonal in the two states which are created by 
the parity odd and even operators 

H-(x) ee H R (x)+H L (x) (2.33) 

H+(x) ee H R (x) - H L (x) . (2.34) 

In fact 

PH±(x) = ±H±(-x) . (2.35) 

The parity odd operator creates the state with lower energy. Thus, it is this the state which 
will have finite energy gap in the continuum limit and which will correspond to the pseu- 
doscalar boson of the continuum Schwinger model. In the continuum one gets this massive 
excitation by applying the vector flux j 1 = j® to the vacuum. The excitation created by H + 
is the scalar excitation which in the continuum limit decouples from the spectrum. Since 
one is interested in the masses of these excitations, one has to compute their energies at zero 
momentum. 

On the lattice the construction of the pseudoscalar excitation at zero momentum is pro- 
vided by 

i N i N i N 

\o> = -^EIW>=-7^E^-W>=-7sE(^ 4 ^ a(+1 +4 + ie- iA ^M >= 

Viv x=1 Viv x=1 Viv x=1 

1 N i N 

J2j 1 (xM>=^=J2j°^M>- (2-36) 



'N^ VN x=l 
One could also consider the scalar massive excitation 

I N I N 

ViV x=1 ViV x=1 



i N i N 

£(4e M V>*+i - iPl+ie- iA A)\ip >= -7= Ei 5 WI^ >l (2-37) 



— VN X=1 
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The states \9 > and \tp > are characterized by a dimer-antidimer configuration; that is by a 
pair of particles followed or preceded by a pair of antiparticles. 

The energy of the state \9 >, at the second order in e, is given by 

E e 2) = E < 0(y)\(H L - H R )A e (H R - H L )\9(x) > (2.38) 

x,y 

where 

A » = if^ (2 ' 39) 

with = E$ + A, A = | and 1 a projection operator onto \9 >. A state obtained 
by applying the hopping Hamiltonian to the state \9 > has a dimer-antidimer pair; as a 
consequence its energy will be raised of 2 A with respect to E^' . [] The computation of E^ 
(see appendix C) leads to 

1 2 

Ef ] = -— N+ — = -AN + 8 (2.40) 

and 

N 2 

< 6\H{K 2 e H h \e >= — - — = 16N - 32 . (2.41) 
Next one should compute 

E t ] = < e{y)\H{K e HlK e H h KeH h \e{x) > + 

-4E < e{y)\HlK 6 H h \9{x) >E< 9{y)\rf h h$H h \0{x) > . (2.42) 
Using the matrix elements given in appendix C, one gets 

Ef ] = -^N - — = 192N - 576 . (2.43) 
The pseudoscalar excitation energy, to the fourth order in e, is then given by 

e 8=(§>+\)+ (~ m + %y + ( i9m - 576 ) e4 • ( 2 - 44 ) 

Using again the matrix element given in appendix C, one easily obtains the energy E v of 
the state \<p >, i.e. 

E v = (— + i) + (-AN + 24)e 2 + (192iV - 1600)e 4 . (2.45) 

Taking the difference between Eg (E^) and one gets the masses of the pseudoscalar 
and scalar excitations , 

p 1 

™ = - + 8e 2 - 576e 4 , (2.46) 
eia 4 



6 One must be careful in applying one more time the hopping Hamiltonian. Since the Coulomb Hamiltonian 
and the hopping Hamiltonian do not commute (see appendix B), there is a state created by appling two times 
the hopping Hamiltonian to the state \6 > which has not a gap of 3A with respect to the ground state, but 
it has a gap of 5A-characterized by three particles followed or preceded by three antiparticles, that is, it is 
created by applying H\ or H\ to the ground state . 



10 



777 I 

4- = ~ + 24e 2 - 1600e 4 . (2.47) 
eia 4 



Eqs.( p.46| , p.47|) provide the behavior of m p and m s for small values of z = e 



2 



4e* a' 



As expected from the continuum theory, the pseudoscalar particle mass is smaller than 
the scalar particle mass. Furthermore, since both quantities are intensive, there is no N 
dependence in ( p.46| , p.47|) . 

Now we compute the dispersion relation of the pseudoscalar and scalar excitations to the 
fourth order in perturbation theory. 

At the zero-th order one has no momentum dependence in the lattice pseudoscalar exci- 
tation since 

E f\p) = ^ E < 0(y)\H \6(x) > e *(-*> = ^ + ^ = ^ + \- ( 2A8 ) 

x,y 

The momentum dependence starts only at the second order in e, where one gets 

E f\p) = < 0{y)\H%H h \d{x) > e^-y^ . (2.49) 

Using the matrix elements given in appendix C, one has 

1 2 

Ef\p) = - — (JV-4) - — cos pa = -AN + 16 - 8 cos pa . (2.50) 
The dispersion relation, at the fourth order in perturbation theory, reads 

E f\p) = < e{y)\H%H%H h A e H h \e{x) > e*-»> + 

x,y 



1 

TV 2 



J2 < e{y)\H%H h \e{x) > e^'rt £ < 6(y)\HtA 2 e H h \9(x) > e*<>*-v) (2.51) 

x,y x,y 

and leads to 

E< i?\p) = ^ N ~ ^ cos2 P a + ^(2cospa-cos 2 pa) - = 

= 192A — 32 cos 2pa + 256(2 cos pa — cos 2 pa) — 800 . (2.52) 

Using ( [2.48] ), Q2.49 ) and ( 2.52 ), the dispersion relation of the lattice pseudoscalar excitation 
is given by 



E e {p) = Ef\p) + e 2 Ef\p) + e*E ( f\p) + ... 



2 " )^ + [^(l - cos2pa) - ^ 

E e (0) + 8(1 - cospa)e 2 + [32(1 - cos2pa) - 256(1 - cospa) 2 ]e 4 . (2.53) 



E o(fl) + ^i 1 ~ cospa)e 2 + [-^i 1 ~ cos2pa) - — (1 - cospa) 2 ]e 4 



Performing a small momentum expansion, one gets 

E e (p) = E e (0) + (4e 2 + 64e 4 )(pa) 2 + . . . (2.54) 
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The same procedure yields 

E v {p) = E^(p) + e 2 E^(p) + e i E^(p) + ...= 

2 1 16 4 

= E e (0) - ^(1 - cospa)e 2 + [— (1 - cos 2pa) + — - — (1 + cospa) 2 ]e 4 = 

= E e (0) - 8(1 - cospa)e 2 + [32(1 - cos2pa) + 1024 - 256(1 + cospa) 2 ]e 4 . (2.55) 

which, for small momentum gives 

E v (p) = E v (0) + (-4e 2 + 72e 4 )(pa) 2 + ... (2.56) 

The p 2 terms are consistent with the bosonic nature of the excitations. 

Eq. (|2.54|) and Eq. (|2.56|) give us the masses of the pseudoscalar and scalar excitations from 



the curvature of the energy momentum relation 



' (4e 2 + 64e 4 ) (2.57) 



2mf e 2 a 3 

= (-4e 2 + 72e 4 ) (2.58 



1 



2m^e 2 a 3 

3 The chiral condensate 

In the continuum Schwinger model, the phenomenon of dynamical symmetry breaking of 
the chiral symmetry, is due to the anomaly. The order parameter is the mass operator 
M(x) = %l)(x)ip(x) which acquires a non zero vacuum expectation value, giving rise to the 
chiral condensate || 

y c =< ib(x)ib{x) >= m c = =. (3.1) 

In this section we compute the lattice chiral condensate to the fourth order in perturbation 
theory. In the staggered fermion formalism (having put particles on even sites and antiparticles 
on odd sites), one has 

^(x)^(x) — - h (3.2) 

a 2 

The lattice chiral condensate may be obtained by considering the mass operator 

1 N 

M = ~Wa EC- 1 )*^. ( 3 - 3 ) 

where the extra minus sign is there to give the same sign to the lattice and continuum chiral 
condensates (see footnote on page 1), and evaluating the expectation value of Q3.3|) on the 
perturbed states \p^ > generated by applying Hh to \if) >. One has 

|JV >= \*l> > +ebi > +e 2 bj > (3-4) 

where 

\pl >= ~H h \il> > (3.5) 

\p% >= ^H h H h \^ > . (3.6) 
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The lattice chiral condensate is then given by 

= <p^\M\p^ > = < > +6 2 < p\\M\p\ > +e 4 < p\\M\p\ > 

XL " <P^\Pi,> " <V# > +e 2 <p\\p\ > +e 4 <V%\V% > 

One gets the following expressions for the wave functions 

< V# > = 1 (3.8) 
<P 1 M> = ^ (3-9) 

<P 2 M> = N( 2A* ^ (3 ' 10) 

and for the M operator 

<^\M\^> = (3.11) 

<p\\M\p\> = -_L(^-2) (3.12) 

<vl\M\p%> = -^(f- 4 )" 2 -^- 3 )- ( 3 - 13 ) 
To the fourth order in e, the lattice chiral condensate is given by 

4 Lattice versus continuum 

In this section we want to extract some physics from the lattice results we obtained; to do 
this, one should compare the answer of the strong coupling analysis of the lattice theory with 
the exact results of the continuum model. For this purpose one should extrapolate the strong 
coupling expansion derived under the assumption that the parameter z = . \ 4 <C 1 to the 

ie L a 

region in which z ^> 1; this region corresponds to the continuum theory since e\a A — ► 0, 
z — > oo. To make this extrapolation possible, it is customary to make use of Pade 
approximants, which allow to extrapolate a series expansion beyond the convergence radius. 
Strong coupling perturbation theory improved with Pade approximants should be compared 
with the continuum theory. 

As we shall see, the gauge invariant strong coupling expansion we propose, provides a 
very accurate estimate of the observables of the continuum theory, already at the second 
order in powers of z. This strongly suggests that, expanding around the gauge, parity and 
charge conjugation invariant ground state \ip >, leads to a perturbative series converging to 



the continuum theory faster then the one used in [[10], |ll], |16|j . 

We first compute the ratio between the continuum value of the meson mass m c = e c / \pn 
and the lattice coupling constant e^, by equating the lattice chiral condensate ( |3.14j ), to its 
continuum counterpart ( [3.1|) 

- - 32z + 1536z 2 = -— m c . (4.1) 
,2 J 2n 
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Of course, eq.(|4.1D is true only when Pade approximants are used, since - as it stands - the 
l.h.s. holds only for z<l, while the r.h.s. provides the value of the chiral condensate to be 
obtained when z ~ oo. Using the relation 



4z 7 e L 

one gets from ( |4.1|) 



^il (4.2) 



m c Az i27r,l 
"iT = ^^2 



-32z+ 1536z 2 ) (4.3) 



As in ref. [11], due to the factor Zi, the fourth power of (|4.3| ) should be considered in order 
to construct a non-diagonal Pade approximant. Since the strong coupling expansion has 
been carried out up to second order in z, one is allowed to construct only the [0, 1] Pade 
approximant for the polynomial written in ( |4.3| ). One gets 

& 4 = i(-) 4 (4.4) 

K e L } 4t 2l e^ l + 256z 1 ' 

Taking the continuum limit z — ► oo one has 

O 4 = W <«> 



Next we compute the same ratio by equating the pseudoscalar mass gap given in fl22HD to its 
continuum counterpart m c 

e\a + 8z - 576z 2 ^j = m c . (4.6) 

Again eq. ([4.6|) is true only when Pade approximants are used. 

Dividing both sides of eq. (|4.6|) by and taking into account that 



t 2 



i 



e L a ={—)*, (4.7) 
one gets 

— = (t- 1 t + 8^- 576z 2 , 4.8 

e L Az 1 v 4 y v 7 

Taking the fourth power, as we did for the chiral condensate equation, and constructing the 
[1,0] Pade approximant for the r.h.s. of eq. (|4.8|) one gets 

(^) 4 = -(— + -) (4.9) 

e L J 4^ 256 T v ! 



One may now take the limit z — > oo, obtaining 



,m c 4 t 2 



"e L ' 8 

Equating eq.(|4.5|) and eq.( [4.10| ) one gets an equation for t which gives 



(4.10) 



t = -J— = 1.049 (4.11) 
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which lies 4.9% above the exact value. It is conforting to see that the lattice theory gives a 
light velocity very close to (but greater than ) 1. Putting this value oft in eq.( |4.5D or eq. ([4. 10|) 
one has 

771 

— = 0.609 (4.12) 

e L 

which lies 7.9% above the exact value -7=. 

It is worth to stress that we can reproduce very well the continuum results, even if we use 
just first order (in z) results of the strong coupling perturbation theory. 

A second way [ffOl to extrapolate the polinomial in (f4.3|) is to equate it to the ratio 



1 + ^ r = 1 - 64z + 3072z 2 (4.13) 



(1 + az) x+ * 

so that one can take the limit z — > 00 in Q4.3p . Setting t = 1 one gets 

77i c \/2n b x 



e L eT a" • 4 



*H 1 



(4.14) 



It is worth to say that this extrapolation works for a large range of values of the parameter 
x between and 6 giving an answer that lies between 6.9% below and 10% above the exact 
answer. In particular for x = 1 one gets ^ = 0.596 which lies 5.6% above the exact value. 
For x = f one gets — = 0.593 wich differs from the continuum value of the 5.1%. For x = 
one has only 0.7% of error. In the same way one can extrapolate the polinomial in 



'1 + bz) x+ i 



1 + az) x 



1 + 32z - 2304z 2 (4.15) 



so that one can then take the limit z — > 00 in (|4.8| ). Setting t = 1 one gets 



- = i (bT + * \ (4-16) 
e L 43 a x 



for x = 1 one gets ^ = 0.617 which lies 9% above the exact value, against the 30% error of 
ref.[|l(|. For x = 0.5 one has ^ = 0.595 which lies only 5.4% above the exact answer. 

In the continuum Schwinger model the ratio between the masses of the scalar and the 
pseudoscalar particles, m s /m p , equals 2, since m s is the lowest eigenvalue of the continuum 
mass spectrum, which starts at 2m p fll |. We shall now evaluate this ratio using the lattice 



expression of m p and m s (p.46|j2.47| ). m s /m p is given by 

m s \ + 24^ - 1600^ 2 



m p 7 + 8z-576z 2 ' ^ 4 ' 17 ^ 



4 



Expanding the r.h.s. of ( [4.17] ) in power series of z 



771 

— = 1 + 64z - 61442 2 , (4.18) 
?7r 



one may perform the [1, 1] Pade approximant 

77l S 1 + 160Z 

m p 1 + 96z 



(4.19) 
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For z — > oo 

— = 1.67 (4.20) 

which lies 16% below the continuum value. Our result coincides with the one obtained in [] 1 
at the same order in the perturbative expansion. 

Taking the [1,1] Pade approximants of eq. ( |2.57| ) and eq. ( |2.58| ) one can compute in an 
independent way the ratio between the masses of the scalar and pseudoscalar particles 



m 



f l + 18z 



(4.21) 



mf 1 — IQz 
For z — > oo 

c 

TYX 

— ^ = 1.125 (4.22) 



mi 



which lies 43% below the exact value. This means that one should go to the next order in 
perturbation theory in eq.( |2.57| ) and eq.( |2.58 ) to reproduce better the continuum limit, since 
in eq.( |2.57| ) and eq.( |2.58| ) one has no zeroth order term. Another test (also suggested in 



TTJ) to check the validity of the lattice computations, may be performed by computing the 



quantity 

D = ~p~r (* 1/4 ™ P ) = 7 + ^4- ■ 

jtt az K ' 4 m/ dz 

This should equal 1/4 if the lattice theory has to reproduce the continuum result. 

Using ( |2.46|) one gets 

s _ 1 1 x,? 7 

D = z , ~ 32z(l - 176z) . (4.24) 

1/4 + 82 v ; v ; 

Constructing the [0, 1] Pade approximant and taking the z — > 00 limit, one has 

D = z 32 z -^= — = 0.182 (4.25) 
l + 176z 11 V ; 



which lies 27% below the desired 0.25. This coincides with the result obtained in [11], using a 
strong coupling expansion up to the order z 4 . The agreement of our results with the continuum 
theory is very encouraging, since only terms up to the order z 2 have been used in the gauge 
invariant strong coupling expansion proposed in this paper. 



5 Concluding remarks 

Our paper is aimed at constructing an improved strong coupling expansion for the lattice 
Schwinger model. We show that, gauge invariance together with the discrete symmetries of 
the model, force the ground state to be a Neel state. Moreover gauge invariance requires the 
ground state to have electric fields related to the charge density operator by the Gauss's law. 
Thus, for large e|, the ground state energy is of order e| rather than ^- |TU| , [TTJ; as a conse- 
quence, the strong coupling expansion we construct provides a very accurate extrapolation to 
the continuum theory already at the second order in z = 4e i a4 . Furthermore, the Coulomb 
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energy of elementary excitations is also affected, since — besides the Coulomb self energy — 
there is also an interaction energy with the charge of the ground state. 

Strong coupling QED with Nl flavors of staggered fermions has been shown to be equiv- 



alent — in any space dimension D — to a quantum SU(Nl) antiferromagnet fll2 |. For the 
strongly coupled lattice Schwinger model, one has equivalence with a one-dimensional an- 
tiferromagnetic Ising spin chain with a Coulomb long range interaction, which admits the 
Neel state as an exact ground state. Due to this, in our approach, chiral symmetry is broken 
spontaneously in the lattice model already at zero-th order in z. We expect that chiral sym- 
metry breaking should persist for all coupling, since the critical coupling for D = 1 should 
be at e L = 0. We conjecture that, for small e 2 , this behavior is a manifestation of a Peierls 
instability — the tendency of a one dimensional Fermi gas to form a gap at the Fermi surface. 
This happens with any infinitesimal interaction. 

In the continuum theory the Schwinger model exhibits the Nambu-Goldstone phenomenon: 
global chiral symmetry is spontaneously broken, but no Goldstone boson appears in the 
spectrum since the local current is anomalous ||. On the lattice there is no anomaly due to 
the Nielsen- Ninomiya theorem 0]. This is obviously true in the staggered fermion formalism 
since 

= E(V^*+i + ipl+M (5-i) 

X 

and 

Q 5 = i[H 3 ,Qs] = , (5.2) 

which leads to no net-production of particles with a given handedness. 

Our analysis shows that in the strong coupling limit, chiral symmetry is spontaneously 
broken at all orders, since the states \ip > and \x > mix only at a perturbative order compara- 
ble with the volume of the system. Being chiral symmetry replaced by the discrete symmetry 
representing invariance under tranlations by one lattice site, effects related to the breaking of 
the chiral symmetry on the lattice should come from the coupling between the two sublattices. 
This is manifest in our strong coupling calculations. 

It will be interesting to extend our gauge invariant strong coupling analysis to include at 
least the case where Nl = 2; in this case, it is expected |T2| that the effective Hamiltonian 
for the model is the antiferromagnetic Heisenberg Hamiltonian. This analysis could provide a 
framework to analyse the properties of the spinon and holon excitations of quantum condensed 
matter models in the context of gauge theories. 



6 Appendix A 

The doubling problem of lattice fermions and a solution to cure it, are reviewed. Let us 
consider the massless Dirac equation in 1 + 1 dimensions in the continuum 



iip = —iad x ip = —i^^d x ip (6.1) 



with ip = I ^ j and a = 75 = u\. Choosing for the plane waves 



V>±0M) = e-* kx - Et h± (6.2) 
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with 

75X± = ±X± (6-3) 

the dispersion relation is 

E{k) = ±k , -oo < k < oo (6.4) 
and the excitations are left- and right-handed particles and antiparticles. 



1> 

each site of a spatial lattice and change the partial derivative d x with a finite difference 



Let us consider the naively latticized Dirac equation. We place the spinor if) = ( ^ j on 



iiP(x) = -— 7 5 [^(x + 1) - ij)(x - 1)] (6.5) 
The dispersion relation is now 

sin(A;a) 

E(k) = ± — — . 6.6 
a 

This is a relativistic dispersion relation for k <C 1 

E(k) ~±k + 0(k 3 a 2 ) (6.7) 
and for ka = it — k'a, with k'a <C 1, 

E(k) ~ Tk' + 0{k' 3 a 2 ) (6.8) 

Let us now turn to the doubling problem of lattice fermions. There are two two-component 
Dirac particles in the continuum limit, and the total chiral charge of the fermions is zero. 
Consider, for example, the excitations associated with the field ip+ on the lattice. There are 
right-movers (k ~ 0) and left-movers (k ~ 7r), and since chirality (helicity for particles) is just 
velocity in 1 + 1 dimensions this lattice field describes a pair of fermions with net chirality 
zero. The doubling problem of lattice fermions is related to the continuous chiral symmetry 
0. In the literature one can find three methods to cure this problem: Wilson fermions, 
SLAC derivatives and staggered fermions. The Wilson fermion formalism lifts the energy at 
the edge of the Brillouin zone. The SLAC formalism uses a nonlocal lattice derivative. In 
what follows we shall review only the staggered fermion formalism that we have adopted in 
this paper. In this formalism one eliminates the unwanted fermion modes by reducing the 
Brillouin zone, i.e. by doubling the effective lattice spacing. In the case of a ci-dimensional 
lattice, one subdivides it into elementary <i-dimensional hypercubes of unit length. At each 
site within a given hypercube one places a different degree of freedom and repeat this structure 
periodically throughout the lattice. Since there are 2 d sites within a hypercube, but only 2a 
components of a Dirac field, we need 2 2 different Dirac fields to reduce the Brillouin zone by 
a factor of |. The lattice Schwinger model in the Hamiltonian formalism |§ with one species 
of staggered fermions ||J] represents the very particular example in which the number of 
sites of the "hypercube" is equal to the number of components of the Dirac field. 

The main idea is to reduce the number of degrees of freedom by using a single component 
Fermi field ip x on each site of the lattice 

{ip x ,ip y } = , {V4 VV> = <W ( 6 - 9 ) 
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so that the lattice Dirac Hamiltonian reads 



H d = ~ E(^+i^ - (6.10) 
and the equations of motion are 

ip x = -i[H D ,if} x ] = T^bPx+i - V'x-i] (6.11) 

If one decomposes the lattice into even and odd sublattices (characterized by x even and odd), 
one can identify a single two component Dirac spinor by associating the upper component 
with even sites and the lower component with odd sites (or vice versa). In this way the 
equations of motion become 

Mx) = -^[il> 2 {x + l)-Mx-l)] (6-12) 

Mx) = -^bl>i(x + i)-Mx-i)]- (6-13) 



When one compares eq.( |6.12| )and eq. ( |6.13| ) with the the continuum Dirac equation, it is easy 



to understand how the staggered fermion formalism works: it avoids the species doubling 
problem by doubling the lattice spacing and so halving the size of the Brillouin zone. 

Let us now discuss the symmetries of the lattice Dirac Hamiltonian ( |6.10| ). It is invariant 
under translations of the spatial lattice by an even number of sites. Translation by two lattice 
spacings is the ordinary continuum translation, whose generator is 



V 



-i J dxip^dxip = -i J dx(ip\d x ipi + i^d^) (6.14) 



and it does not mix upper and lower spinor components. Consequently, one expects the lattice 
generator not to mix the two sublattices 

p = -i ^2(ipl +2 ipx + i/>U>x+2) (6-15) 



The Hamiltonian (|6.10|) is also invariant under translations by an odd number of sites. 
The lattice generator of this symmetry is 

^ = £(^+1+^+1^) (6-16) 

X 

which in the continuum would be T5 = / dxip'^ip. The matrix 75 applied to the bispinor 
interchanges upper and lower component, this on the lattice corresponds to translation by 
one lattice site. Therefore, the staggered fermion formalism has a discrete chiral symmetry, 
corresponding to chiral a rotation of | radiants. 



Let us consider now the naive continuum limit of the lattice Dirac Hamiltonian ( |6.10|) . 
We consider an elementary "plaquette" of the lattice as the segment of length a with two 
sites generated by taking a site of even coordinates and adding to it the lattice spacing a. We 
also decompose the lattice into two sublattices generated by taking a site of the elementary 
segment and translating it by all even multiples of the lattice spacing. We label the fermions 
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which reside on the sublattice of each of the site of the elementary segment as ip a (x) with 
a = if one is on the first site and a — 1 if one is on the second site. In momentum space 
the Hamiltonian is 

H D = [ dk^ a (k)T l a(3 sin ka^^k) (6.17) 

where Qb = {k : — ^ < k < ^} is the Brillouin zone of the even sublattice. The momentum 
space fermions have the anticommutator 

{Mh),il>l(k , )} = 6 a ( j 6 k , k , (6.18) 

and the Dirac tensor is 

= alp (6.19) 

The spectrum of the Dirac operator is 

sin fca . . 

E(k) = ± (6.20) 

a 

and one notes that, to set the Fermi level of the fermions at the degeneracy point where the 
two branches of the spectrum meet, it is necessary that the fermion states are exactly half 
filled. This is also required for charge-conjugation invariance, or particle- hole symmetry of 
the vacuum state. 



7 Appendix B 

In this appendix the energies of the excitations generated by applying the hopping Hamil- 
tonian Hh to the ground state \ip >, are evaluated. The action of Hh on the state \ip > 
generates a dimer-antidimer pair. All the states with one dimer-antidimer pair, and with two 
dimer-antidimer pairs (obtained applying Hh twice), have the same energy gap with respect 
to the ground state. However, not all the states obtained applying (H h ) n , with n > 3, to \ip > 
have the same energy. There are, for example, two possible energies, corresponding to a gap 
of 3 A = 3e|a/4 and 5 A, for the states with three dimer-antidimer pairs. 

The states with a dimer-antidimer pair are given by 

\S 1 R , L >=H R , L (y)\^>. (7.1) 
Applying on these states the unperturbed Hamiltonian we get 

H u \s R , L > = ^E(^±M 2 I^>= 

A X 

= K + £ |^E(-^^ + M]l^> (7.2) 

where y must be even if we apply H R and must be odd if we apply H L . We thus see that all 
the states ^ > have a mass gap of A. 

The states with two dimer-antidimer pairs are given by 

|S(W,L) >= H R , L (z)H R:L (y)\iP > (7.3) 
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Hu\Sf R , L){RiL) >= ^ J2(E X ± Sxz ± 5 xy ) V > (7.4) 

It is easy to show that all these states have an energy gap of 2 A. 

Based on (|7.2|) and (|7.4|) one might expect that the states with three dimer-antidimer pairs 
have an energy gap of 3 A. This is actually true for all the states of this type except for the 2 
states with three particles followed or preceeded by three antiparticles. These states have an 
energy gap of 5A. If one considers the state 



\Srrr >= H R {w)H R {z)H R {y)\^ >, (7.5) 



one easily obtains that 



6^ Qj 

H u\S RRR > = Yl E * + 5 XW + 8xZ + 5 xy ) 2 \lp > = 

^ X 

e\a (— l) x (— l) x (— l) x 
= i E i> + -|- XX —$xw — S xg — 8 xy + 

^&xw&xz ^&xw&xy ^&xz&xy 

+ S xw + 5 XZ + 5 xy )]\ip > . (7.6) 

To generate a state with three antiparticles followed by three particles, one must have 

z = V - 2 (7.7) 
w = y-l (7.8) 

The energy gap for this state is 5 A. The same result is obtained for the state \S\ LL >= 
Hl{w)Hl{z)Hl{ii)\iP > with z = y + 3 and w = y + 2, which contains three particles followed 
by three antiparticles. 

The reason for these different behaviours, lies on the fact that the commutation relation 
between H and is 

[Ho, H h ] = ^H h + %Y,~ T,M + ie lA ^k ~ 4e~ tA Mp 3 (7.9) 

j>k k>j 

If the second term in the r.h.s. was absent, all the states with n dimer-antidimer pairs would 
have an energy gap of nA. 



8 Appendix C 

In this appendix we provide the matrix elements needed in the strong coupling expansion 
expression of the mass gap. The six matrix elements that arise in the computation of the 
second order bosonic excitation energies, are 

M 1 = < ^\H R (y)H R A H L H L (x)\4j >= 

= ~(f -1)^-1(1-^) (8.1) 
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M 2 = <il>\H L (y)H L A e H R H R (x)\i/;>-. 



Ma 



M R 



IN 1 

< ^\H R (y)H L A e H R H L (x)\^ >= 
l.N 1 

< ^\H L (y)H R A e H L H R (x)\tlj >= 
IN 1 

- < ^\H R (y)H L A e H L H R (x)\^ >-- 

- < ijj\H L (y)H R A H R H L {x)\i) >- 



(8.2) 
(8.3) 
(8.4) 
(8.5) 
(8.6) 



where Mi and M 3 are defined for x and ?/ odd, M 2 and M 4 for x and y even, M 5 for x even 
and ?/ odd, and M 6 for x odd and y even. 

To compute the fourth order energy gap one needs also the matrix elements M-, i = 1, 6 
obtained by replacing Aq with A@ in Mj 

1 JV 1 

M[ = M' 2 = M'^ = M' 4 = —{- - l)5 xy + — (1 - 5 xy ) (8.7) 

M 5 = M 6 = - — (£j,,x-i + <Wi) - ^2 (l-^-i -<W0 (8-8) 

with the same constraints on x and y as before. 

One must then evaluate the twenty matrix elements arising when three H R and three Hi 
are combined in all possible ways in the expression in the first line of eq. ( |2.42|) 

< 6{y)\{H L - H R )A e (H L - H R )A e (H R - H L )A e (H R - H L )\6(x) > (8.9) 

Since hermiticity of the Hamiltonian requires that the elements obtained changing H R with 
Hl are equal one has only to compute the ten matrix elements, 



M'l = 


— y M'l 

N ^ 1 

x ,y 


x,y) = 


N ^ 

x,y 


< ^\H R (y)H R A e H R A e H L A e H L H L 


[x)\ip > 


(8.10) 


M'i = 


x,y 


>,y) = 


N ^ 

x >y 


< ^\H R (y)H L A e H L A e H R A e H R H L 


[x)\ip > 


(8.11) 


M'l = 




x,y) = 


N ^ 

x >y 


< 4>\H R (y)H R A e H L A e H R A H L H L 


[x)\ip > 


(8.12) 


M'l = 
M'l = 


-Ym'I 

N 4 

x ,y 

x,y 


x,y) = 
x,y) = 


N ^ 

x >y 

N ^ 

x ,y 


< i>\H R (y)H R A e H L A e H L A e H R H L 

< ^\H R {y)H L A e H R A e H R A e H L H L 


[x)\ij) > 
[x)\ip > 


(8.13) 
(8.14) 


M'l = 


x >y 


x,y) = 


N ^ 

x,y 


< ^\H R (y)H L A e H R A e H L A e H R H L 


[x)\ip > 


(8.15) 
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M!} = 


1 

N 


x,y 


-^E < mR(y)H R A 9 H L AeH L A e H L H R (xM > 

x,y 


M'i = 


1 

N 




-^E < mR(y)H L A e H R AeH L A e H L H R (x)\ij > 

x,y 


Ml = 


1 

TV 


x,y 


~^E < mR{y)H L AeH L AeH R A e H L H R {x)\i, > 

x,y 


<o = 


1 

N 


£M(' (x,y) = 

x,y 


-^E < mR{y)HLAeH L A e H L A e H R H R {x)\i) > 

x,y 



where M"(x, y), M'l (x, y) are defined for x and y odd, M" (x, y), Mf (x, y) are defined 
for x even and y odd. The results are 

K(x,y) = -^6 XtV -± i (N-2)(N-4)6 x , v + 

- ^fc + 2 + ^2)-^(iV-4)(l-^) (8.20) 

M?(x,y) = -JL(N-4)(N-6)5 X „ (8.21) 

M 3 '(x,y) = M':(x,y) = Ml(x,y) = Ml(x,y) = -^(N-A)(N-Q)5 x , y + 

M?(x,y) = K(x,y) = M^(x,y) = Mf (x,y) = 

= 2^( JV - 6 )( 1 "^-i- W) ( 8 -23) 

and 

M i = "^-^( iV - 2 )( iV - 4 ) ( 8 - 24 ) 

M £ = A# = M£ = M£ = M£ = -^(JV-4)(iV-6) (8.25) 

M ; = M» = MS = = +±(N - 4)(N - 6) (8.26) 
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